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16 Abstract 

17 Snow cover over the Northern Hemisphere plays a crucial role in the Earth’s hydrology 

18 and surface energy balance, and modulates feedbacks that control variations of global 

19 climate. While many of these variations are associated with exchanges of energy and 

20 mass between the land surface and the atmosphere, other expected changes are likely to 

21 propagate downstream and affect oceanic processes in coastal zones. For example, a large 

22 component of the freshwater flux into the Arctic Ocean comes from snow melt. The 

23 timing and magnitude of this flux affects biological and thermodynamic processes in the 
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1 Arctic Ocean, and potentially across the globe through their impact on North Atlantic 

2 Deep Water formation. 

3 

4 Several recent global remotely sensed products provide information at unprecedented 

5 temporal, spatial, and spectral resolutions. In this article we review the theoretical 

6 underpinnings and characteristics of three key products. We also demonstrate the 

7 seasonal and spatial patterns of agreement and disagreement amongst them, and discuss 

8 current and future directions in their application and development. Though there is 

9 general agreement amongst these products, there can be disagreement over certain 

10 geographic regions and under conditions of ephemeral, patchy and melting snow. 

11 
12 
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1 1. Introduction 

2 

3 Snow covers a considerable portion of Northern Hemisphere lands during winter. It is the 

4 component of the cryosphere with the largest seasonal variation in spatial extent. In fact 

5 accumulation and rapid melt are two of the most dramatic seasonal environmental 

6 changes of any kind on the Earth’s surface (Gutzler and Rosen 1992, Robinson and Frei 

7 2000, Robinson et al. 1993). In the Southern Hemisphere, outside of Antarctica and its 

8 surrounding ice shelves and sea ice, snow is generally limited to smaller regions such as 

9 the Andes, Patagonia and the southern Alps of New Zealand (Foster et al. 2008). On 

10 decadal time scales, snow variations over Northern Hemisphere lands have also been 

11 considerable (Barry et al. 1995, Brown 2000, Brown and Braaten 1998, Derksen et al. 

12 2004, Frei et al. 1999, Mote 2006, Mote et al. 2005, Ye et al. 1998), with declines in 

13 spring associated with warmer conditions (Brown et al. 2010, Groisman et al. 1994, IPCC 

14 2007, Feathers and Robinson 1993). Recent reports on changes in the Arctic environment 

15 cite snow as one of the critical variables (ACIA 2004, AMAP 2011). The expectation 

16 during the 21 st century is that changes will be increasingly dramatic (Frei and Gong 2005, 

17 Raisanen 2007, Ye and Mather 1997) and spatially and temporally complex (Brown and 

18 Mote 2009, Nolin and Daly 2006). 

19 

20 While large scale changes in snow cover are useful as indicators of climatic variations, 

21 snow also affects other components of the Earth system at a variety of scales. By virtue 

22 of its radiative and thermal properties which modulate transfers of energy and mass at the 

23 surface-atmosphere interface (Zhang 2005), snow affects the overlying atmosphere 


Page 3 of 52 



1 (Barry 2002, Barry et al. 2007, Cohen 1994, Ellis and Leathers 1999, Mote 2008, Walsh 

2 1984) and thereby plays an important role in the complex web of feedbacks that control 

3 local to global climate. For example, because of the high albedo of snow, its presence can 

4 change the surface energy balance over land and ice and therefore affect climate (i.e. the 

5 snow-albedo feedback). Snow also modulates the hydrologic cycle (Dyer 2008, Graybeal 

6 and Leathers 2006, Leathers et al. 1998, Todhunter 2001); influences ecosystem 

7 functioning (Jones et al. 2001); and is a significant resource for many mid latitude 

8 populations and for populations whose water is derived from mountainous and northerly 

9 cold regions (Barnett et al. 2005, Barry et al. 2007). Snow observations are critical for the 

10 validation of climate models (Foster et al. 1996, Frei et al. 2003, Frei et al. 2005, MacKay 

11 et al. 2006, Roesch et al. 1999). 

12 

13 With regards to the freshwater flux to the ocean, the role of snow is to modulate seasonal 

14 timing, and in some cases the amount, of discharge into the oceans. While this can affect 

15 coastal systems across mid-latitudes, of particular relevance is the fresh water flux into 

16 the Arctic basin. The drainage area into the Arctic Ocean is ~1.5 times the surface area of 

17 the Arctic Ocean itself (Peterson et al. 2002) and river runoff is the largest source of 

18 freshwater input into the Arctic basin (Arnell 2005, Miller and Russell 2000). Much of 

19 Arctic precipitation is derived from snow fall, and much of the river runoff is derived 

20 from snow melt. During the past century, both high latitude precipitation (Zhang et al. 

21 2007) and river runoff to the Arctic basin have increased; both are expected to increase 

22 further in a warming climate (Peterson et al. 2002), although the rates of change and 

23 relative impacts on ocean circulation vary spatially (Rennermalm et al. 2007). 
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1 


2 The studies described above do not include all the possible nonlinear feedbacks in which 

3 snow plays a role in the Arctic environment (Hinzman et al. 2005). For example, due to 

4 the insulating effect of snow cover, changes in the timing of snow onset or disappearance, 

5 or changes in the amount of snow, may influence the state of the underlying permafrost, 

6 which has been warming for decades (Romanovsky et al. 2010) and which is expected to 

7 deteriorate during this century (Lawrence and Slater 2005) and may further increase the 

8 freshwater flux. Thawing permafrost may also result in a significant release of carbon to 

9 the atmosphere as the result of microbial decomposition of currently frozen organic 

10 carbon (Schuur et al. 2008). According to Betts et al. (2000) the expected expansion of 

11 the boreal forest may lead to both negative feedbacks (an additional carbon sink) and 

12 positive feedbacks (an albedo decrease) on global climate, and the net effect will be a 

13 positive feedback with increased warming. The feedbacks between snow, permafrost, and 

14 freshwater flux to the Arctic Ocean associated with these processes are poorly understood 

15 (Francis et al. 2009, Rawlins et al. 2010). 

16 

17 While an increased freshwater flux to the Arctic has potential effects on thermodynamic 

18 and ecological processes in the coastal zone, perhaps most importantly such increases 

19 have been shown in the past to diminish or completely halt the formation of North 

20 Atlantic Deep Water (NADW) (Rahmstorf 2000). This occurs because freshwater export 

21 to the North Atlantic Ocean, the region of NADW formation, decreases surface water 

22 density. Model simulations suggest that the magnitude of expected runoff changes during 

23 this century may approach critical thresholds for NADW formation (Arnell 2005, Miller 
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1 and Russell 2000, Peterson et al. 2002). In a recent study, NADW formation as well as 

2 permafrost degradation and changes to the tundra and boreal forest ecosystems (all of 

3 which can be affected by snow, and all of which can affect the freshwater flux to the 

4 ocean) have been listed among the potentially critical components of the Earth system 

5 that may be in danger of approaching “tipping points” (Lenton et al. 2008). Thus, 

6 accurate monitoring of high latitude snow remains an essential goal. 

7 

8 Because of the large extent of terrestrial snow cover and the difficulties in obtaining 

9 ground measurements over cold regions, remote sensing represents an important tool for 

10 studying snow properties at regional to global scales. In recent years, advances in satellite 

11 capabilities, as well as in algorithm development, have led to improved monitoring of 

12 snow across the globe. The purpose of this article is to review the current generation of 

13 satellite-derived global snow observations that has become available during the first 

14 decade of the twenty first century, with emphasis on land surfaces of the Northern 

15 Hemisphere. Theoretical considerations for the remote sensing of snow, and key products 

16 are discussed. 

17 

18 2. Theoretical Background 

19 

20 Due to the nature of interactions between snow cover and electromagnetic radiation of 

21 different frequencies, snow can be distinguished from other terrestrial surfaces using 

22 satellite observations based on a number of different active and passive techniques 

23 (Dozier 1989, Nolin 2010). The two types of instruments used for monitoring global 
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1 scale snow variations rely on either (1) a combination of the visible and infrared, or (2) 

2 microwave, portions of the electromagnetic spectrum (Hall et al. 2005, Matzler 1994, 

3 Rango et al. 2000, Scherer et al. 2005, Schmugge et al. 2002). These methods are limited 

4 by a number of factors, such as clouds, forest cover fraction, terrain heterogeneity and 

5 precipitation. For example, interpretation of visible and infrared as well as passive 

6 microwave images can be difficult where complex terrain causes considerable spatial 

7 variation within each remotely- sensed footprint of snow depth, surface characteristics, 

8 and satellite viewing angles. Nevertheless, products based on these observations have 

9 been vital for monitoring snow and for our understanding of the role of snow in the Earth 

10 system. Though global active microwave data (e.g., QuikSCAT) can also be used to 

11 study snow extent and depth at relatively large spatial scale (Tedesco and Miller 2007a, 

12 b), data are available only from 1999 through 2009 (when the satellite failed well past its 

13 expected lifetime; see http://www.jpl.nasa. gov/news/news.cfm?release=2009- 175 

14 downloaded November 2011). In contrast, passive microwave data have been available 

15 since the late 1970s, and continue to be available. At regional scales, airborne data can 

16 also be collected before and after the snow falls to study the attenuation introduced by the 

17 snow pack on naturally emitted gamma radiation (Carroll 1987). However, the data 

18 collected with this method have low temporal resolution (seasonal scale) and cannot be 

19 used for global scale studies. Consequently, we focus our analysis on snow parameters 

20 estimated by means of visible and infrared and passive microwave sensors. 

21 

22 2.1 Visible and Near-Infrared 

23 
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1 Snow extent (i.e. presence or absence of snow, regardless of snow amount) is, in many 

2 circumstances, relatively straightforward to observe using visible observations because of 

3 the high albedo of snow (up to ~ 80 % or more in the visible part of the electromagnetic 

4 spectrum) relative to most land surfaces. However, limitations exist. First, visible 

5 imagery is limited to that portion of the surface illuminated by sunlight; thus darkness 

6 and low illumination are problematic. Second, clouds impede visible evaluation in two 

7 ways. All but the thinnest clouds reflect a significant portion of visible radiation, 

8 preventing any visible radiative information about the surface from reaching the satellite. 

9 And, because the albedos of clouds and snow are often similar, the discrimination 

10 between cloud-covered and snow covered surfaces can be difficult. However, near- 

11 infrared bands can be used to distinguish between snow and most clouds because the 

12 near-infrared reflectance of most clouds is high while the near-infrared reflectance of 

13 snow is low. 

14 

15 Third, vegetation can obstruct visible and infrared information about snow from reaching 

16 the satellite sensor. Forest canopies protrude above the snow pack, lowering the surface 

17 albedo (Robinson and Kukla 1985) and partially or completely obscuring the underlying 

18 surface, making it difficult to determine snow extent or amount (Chang et al. 1996, 

19 Derksen 2008, Goita et al. 2003, Klein et al. 1998, Nolin 2004). 

20 

21 Lastly, surface heterogeneity can play a role the interpretation of visible and infrared 

22 imagery in a number of ways. Of particular relevance to the monitoring of high latitude 

23 snow is the presence of numerous frozen lakes in Arctic regions, which may contribute to 
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1 the overestimation of snow covered area from visible and infrared based imagery during 

2 periods when lakes remain frozen after the snow has melted on adjacent land surfaces 

3 (Derksen et al. 2005a, Frei and Lee 2010), at least when high resolution land surface data 

4 sets are not used to filter out the signal from lake surfaces. Passive microwave based 

5 estimates of SWE may be underestimated due to the presence of lakes (Derksen et al. 

6 2005a, Rees et al. 2006). On the other hand, surface heterogeneity may assist in the 

7 interpretation of snow-covered versus snow-free ground, and of snow-covered versus 

8 cloud-covered scenes, when trained analysts are mapping snow extent using visible 

9 imagery. 

10 

11 2.2 Passive microwave 

12 

13 Because snow grain dimensions can be similar to microwave wavelengths, snow is 

14 efficient at scattering the microwave radiation naturally emitted from the Earth’s surface 

15 (Matzler 1994). Therefore, microwave emission from a snow covered surface is 

16 diminished relative to a snow-free surface, and the presence of snow can frequently be 

17 identified (Chang et al. 1976, Grody 2008, Hall et al. 2005, Matzler 1994, Tait 1998, 

18 Tedesco and Kim 2006). Furthermore, because under ideal circumstances the amount of 

19 scattering is proportional to the number of snow grains, microwave instruments offer the 

20 possibility of estimating the mass per unit area of water in the snow pack, which is often 

21 measured as snow water equivalent (SWE). In contrast to visible and infrared, passive 

22 microwave does not depend on the presence of sunlight and thus provides an alternative 

23 at high latitudes; and, passive microwave is largely (but not completely) transmitted 
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1 through non-precipitating clouds, offering the potential to estimate snow cover under 

2 many cloudy conditions that preclude visible and infrared observations. In practice, 

3 research using passive microwave exploits the fact that microwave scattering by ice 

4 crystals is frequency-dependent: higher frequencies within the microwave portion of the 

5 spectrum are scattered more efficiently than lower frequencies, enabling the use of two or 

6 more frequency bands to estimate SWE (Chang et al. 1987, Derksen 2008, Derksen et al. 

7 2005b, Grody and Basist 1996). Other methods have also been evaluated such as one 

8 based on the inversion of a snow emission model (e.g., Pulliainen and Hallikainen 2001). 

9 Clifford (2010) provides a review of global estimates of snow water equivalent from 

10 passive microwave. 

11 

12 Limitations to the monitoring of snow using passive microwave sensors are due to a 

13 variety of factors. One major limitation is the presence of liquid water in the snow pack, 

14 the microwave emission from which masks the snow signal and inhibits the ability of 

15 microwave sensors to detect wet snow. Also, because of the relatively weak microwave 

16 signal emitted by terrestrial surfaces, microwave sensor footprints are necessarily large 

17 (~25 km). Uncertainties in snow depth and SWE estimates are associated with the 

18 physical structure of snow packs (ice lenses, grain size variations and vertical 

19 heterogeneity) which vary in space (Chang et al. 1976, Sturm et al. 1995) and time 

20 (Langham 1981) and can alter the scattering and emission characteristics of the snow 

21 pack. Snow pack metamorphosis, which in the Arctic region typically results in a layer of 

22 depth hoar (with large crystal size) near the bottom of the pack, results in more efficient 

23 microwave scattering. Thus, a signal change measured at the satellite sensor due to snow 
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1 metamorphosis can mimic a signal change due to a change in SWE. Vegetation in and 

2 above the snow pack emits microwave radiation, and can confound any detection 

3 algorithm (Chang et al. 1996, Foster et al. 1997, Tedesco et al. 2005). 

4 

5 Finally, as a snow pack reaches a certain critical depth the relationship between snow- 

6 amount and MW brightness temperature reverses (Derksen 2008, Markus et al. 2006, 

7 Matzler 1994, Schanda et al. 1983). When SWE exceeds -150 mm emission by the snow 

8 pack of microwave band radiation is greater than scattering, resulting in a positive 

9 relationship between SWE and brightness temperature. This is an additional source of 

10 uncertainty in SWE retrievals for deep snow packs. 

11 

12 2.3 Remote sensing of snow in complex terrain 

13 Some of the difficulties inherent in the interpretation of remotely sensed images are 

14 exacerbated in regions with complex terrain (Dozier 1989). Due to variability of slope, 

15 aspect, and land cover, the local solar illumination angle varies within one satellite 

16 footprint. In fact, due to co -registration differences between an image and a digital 

17 elevation model, illumination angles, and therefore reflectance characteristics, are often 

18 unknown. To address such issues. Painter et al. (2009) developed the MODSCAG model, 

19 which estimates mean grain size and fractional snow cover from MODIS data using 

20 linear spectral mixture analysis and a library of reflectance characteristics of different 

21 surface types. This model has relatively small errors, and could potentially be applied 

22 globally, but so far has been validated mostly in regions of complex terrain. 

23 
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1 A recent study of different algorithms for estimating SWE from passive microwave 

2 radiances in a basin with complex terrain in the Canadian Rockies finds that the 

3 traditional algorithms which are based on brightness temperature differences in difference 

4 wavelength intervals (as discussed above) are less accurate than Artificial Neural 

5 Network (ANN) techniques which can be trained on observations from surface stations 

6 (Tong et al. 2010a). Unfortunately, due to the limited distribution of stations for training 

7 the ANN in their test region, they are unable to accurately estimate spatial variations of 

8 SWE. Tong and Velicogna (2010) and Tong et al. (2010b), using surface station 

9 observations and MODIS imagery across the Mackenzie River Basin, determine that the 

10 minimum or threshold SWE value estimated from passive microwave observations that 

11 can be used to determine the presence / absence of snow varies from sub-basin to sub- 

12 basin, depending on topography and vegetation cover. Nevertheless they find useful 

13 information in remotely sensed SWE values for hydrologic monitoring. As these studies 

14 indicate, the estimation of snow characteristics in complex terrain from remotely sensed 

15 imagery is an important and cutting edge field of study. At this time, these techniques 

16 have not been incorporated into global products, and are not addressed further in this 

17 paper. 

18 

19 2.4 Comparison and evaluation of products 

20 

21 When two products disagree, which is “correct?” Is either one of them “correct?” Two 

22 key impediments to a conclusive evaluation are that there is no perfect “ground truth,” 

23 and that the answer depends on spatial scale. The most obvious method of testing the 
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1 veracity of remotely sensed (or other gridded) products is by comparison to surface 

2 reference observations. However, there exists considerable contrast between surface, or 

3 in situ, and remote snow observations with regards to the snow pack properties that can 

4 be measured, their spatial and temporal resolutions and domains, and the methods 

5 employed to make measurements (Brown and Armstrong 2008, Goodison et al. 1981). 

6 

7 Even in regions with surface observations, validation may be difficult because of the 

8 contrasting spatial scales of surface and remotely sensed observations (Brubaker et al. 

9 2005, Chang et al. 2005). Brubaker et al. (2005) discuss the difficulties in comparing 

10 point measurements to spatially integrated satellite retrievals, especially in areas of sparse 

1 1 station networks, which are typically at high elevations and northerly regions (exactly the 

12 areas where snow is most prevalent). They find that there is no single accepted method to 

13 perform validation of remotely sensed snow products. Chang et al. (2005) provide an 

14 informative review of how varying station densities and different satellite footprints are 

15 not equally spatially representative, and how the differences can complicate evaluations 

16 and comparisons of different products. They employ geostatistical techniques, as 

17 suggested by Kelly et al. (2004), to quantitatively define the spatial density of station 

18 observations required to provide sufficient information for validation studies. MODIS has 

19 been found to compare well with station based observations as well as with the National 

20 Operational Hydrologic Remote Sensing Center products (Hall and Riggs 2007), but 

21 Riggs et al. (2005) show that even between different versions of MODIS snow products, 

22 analyses at different spatial resolutions can provide conflicting results in some cases, due 

23 to both the resolution differences and the averaging method. 
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1 


2 Despite the inherent difficulties, comparative studies to date have drawn some useful 

3 conclusions (Armstrong and Brodzik 2001, Basist et al. 1996, Bitner et al. 2002, Brown 

4 et al. 2007, Brown et al. 2010, Derksen et al. 2003, Drusch et al. 2004, Foster et al. 1997, 

5 Mialon et al. 2005, Mote et al. 2003, Romanov et al. 2002, Savoie et al. 2007, Tait and 

6 Armstrong 1996). For example, evaluations of NOAA visible and infrared versus passive 

7 microwave products find more disagreement during fall and spring than during mid- 

8 winter, with particular differences under forest canopies, over complex terrain, in areas of 

9 persistent clouds, patchy snow, and wet snow (Armstrong and Brodzik 2001, Basist et al. 

10 1996). Over the Tibetan Plateau these products often disagree year-round (Armstrong and 

11 Brodzik 2001, Savoie et al. 2007). Several recent studies identify differences between 

12 remotely sensed products and surface observations over North America during the spring 

13 ablation season (Brown et al. 2007, Brown et al. 2010, Frei and Lee 2010). 

14 

15 3. Snow Products 

16 

17 A number of digital products that are based on remote observations are available. The 

18 two visible and infrared based suites of products that are most widely used for large-scale 

19 climate research are from: (1) the Interactive Multisensor Snow and Ice Mapping System 

20 (IMS) (section 3.1) and (2) the suite of products derived from the Moderate Resolution 

21 Imaging Spectroradiometer (MODIS) (section 3.2). IMS is the most recent version of a 

22 product that dates back to the 1960s (Matson and Wiesnet 1981). IMS mapping of snow 

23 extent has relied primarily on visible and near infrared imagery, but includes data and 
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1 information from a number of sources. As discussed in more detail below, the key feature 

2 that distinguishes IMS from other products is human involvement in the analysis, which 

3 is required for operational purposes. 

4 

5 The MODIS instrument, which is used to observe a number of geophysical variables 

6 including snow, flies on NASA’s Earth Observing System (EOS) Terra satellite, 

7 launched in 1999. A near-twin MODIS instrument is also flying on board the Aqua 

8 platform, which was launched in 2002 (Aqua and Terra have afternoon and morning 

9 equatorial crossing times, respectively). Aqua also hosted the Advanced Microwave 

10 Scanning Radiometer - Earth Observing System (AMSR-E) until its failure in October 

11 2011. AMSR-E measured the naturally emitted radiation in the microwave region at 5 

12 different frequencies (6.9, 10.7, 18.7, 23.8, 36.5 and 89 GHz) at both vertical and 

13 horizontal polarizations. 

14 

15 The IMS and MODIS snow algorithms both rely primarily on near-polar orbiting 

16 satellites, from which daily images are available at high latitudes. Other algorithms that 

17 have been suggested (Romanov et al. 2003, Siljamo and Hyvarinen 2011) rely on 

18 geostationary satellites, which have the advantage of higher temporal resolution, but have 

19 poor spatial resolution. 

20 

21 3.1 visible and near infrared based products 

22 

23 3.1.1 The Interactive Multisensor Snow and Ice Mapping System ( IMS ) 
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2 The data set that has historically been the most widely used for the operational mapping 

3 and climatological analysis of large-scale snow extent (not depth or water equivalent) 

4 was produced by the US National Oceanic and Atmospheric Administration (NOAA) 

5 National Environmental Satellite and Data Information Service (NESDIS), but has been 

6 transferred to the National Ice Center (NIC), which is jointly supported by NOAA, the 

7 US Navy, and the US Coast Guard. This product has been based primarily on visible and 

8 near infrared observations, and covers the period from late 1998 to present (Ramsay 

9 1998), with the precursor maps beginning in 1966, constituting the longest remotely 

10 sensed environmental time series that has been derived in a near-consistent fashion 

11 (Helfrich et al. 2007, Matson and Wiesnet 1981, Robinson et al. 1993). The term near- 

12 consistent is used because, due to changing operational requirements and evolving 

13 technical capabilities, this product has undergone, and continues to undergo, 

14 improvements and refinements (Helfrich et al. 2007, Ramsay 1998, 2000) as summarized 

15 briefly here. The two reasons for this product’s importance - as operational input into 

16 atmospheric forecast models and as a long-term climatic record - are also discussed. 

17 

18 Although a number of improvements and corrections in the production of the NOAA 

19 product occurred in the earlier years, the biggest methodological change was 

20 implemented in the late 1990s. Until that time, NOAA snow maps were produced on a 

21 weekly basis by trained meteorologists who would visually interpret photographic copies 

22 of visible band imagery, and manually produce maps that would subsequently be 

23 digitized with spatial resolution between 150 km and 200 km. In 1997 NOAA began 
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1 producing snow maps using the IMS, with improved spatial (24 km) and temporal (daily) 

2 resolutions. IMS is operated by trained analysts who produce a daily digital product 

3 utilizing Geographic Information System technology and incorporating a variety of, and 

4 an ongoing expansion of, technological capabilities as well as sources of information. 

5 Since 1999, when weekly manual mapping was discontinued, daily IMS maps have been 

6 produced (Ramsay 1998, Robinson et al. 1999). Technological advancements since 1999 

7 have led to even higher resolution (4 km) snow mapping (Helfrich et al. 2007). 

8 

9 IMS produces estimates of snow extent across the globe every day, regardless of the 

10 presence of clouds. This is possible primarily for two reasons. First, analysts use sources 

11 of information other than visible and near infrared imagery. Second, because IMS 

12 analysts can loop through sequential images, their ability to evaluate scenes is based on 

13 an integration of information from both spatial and temporal perspectives. Thus, a key 

14 feature of the IMS product is that human judgment as to which data sources are most 

15 reliable in different conditions and regions, and as to the final evaluation of where the 

16 snow is, remains an integral part of the process, and one of the strengths of the IMS 

17 product. IMS also includes sea ice extent, which is not discussed in this report. Figure 1 

18 shows an example of a daily IMS snow map in its original projection. 

19 

20 It is difficult to optimize this product for both of its two main uses. Its primary purpose is 

21 to provide input to atmospheric forecast models. For this purpose, continued product 

22 improvements are advantageous. As a record for evaluating long term environmental 

23 change, however, the value of any product is diminished if methodological changes 
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(including those that provide more accurate maps) result in temporal inconsistencies in 
the data set that might be difficult to distinguish from actual variations in snow extent. To 
maintain product continuity and a viable long-term record, IMS continues to produce a 
coarse (24 km) resolution version of the data set. And, in collaboration with NOAA, the 
Rutgers University Global Snow Lab is producing a climate data record in which 
inconsistencies between the earlier maps and the IMS product (in addition to 
inconsistencies during the weekly map era) are accounted for, and can therefore be used 
for ongoing analyses of historical variations (Robinson and Estilow 2008). 

3.1.2 The Moderate Resolution Imaging Spectroradiometer (MODIS) 

The MODIS sensor measures radiation in 36 spectral bands, including the visible, near 
infrared, and infrared parts of the electromagnetic spectrum. The suite of MODIS snow 
cover products, available since 2000, are derived using a fully-automated algorithm that 
provides good spatial resolution (500m), cloud detection, and frequent coverage (daily at 
mid to high latitudes) (Hall and Riggs 2007, Hall et al. 1995, Hall et al. 2002, Riggs et al. 
2006). The MODIS snow-mapping algorithm uses a normalized difference between 
MODIS band 4 (5.45-5.65 mm) and 6 (1.628-1.652 pm) and many additional spectral 
and threshold tests. In forested areas the threshold is changed based on results of a 
canopy reflectance model, using both the Normalized Difference Vegetation Index 
(NDVI) and Normalized Difference Snow Index (NDSI) (Klein et al. 1998). A thermal 
mask is also included to remove erroneous “snow” over locations where the presence of 
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snow is considered to be implausible. See Riggs et al. (2006) and Riggs and Hall (in 
press) for a description of the algorithm and recent upgrades. 

NASA provides a hierarchy of snow products based on MODIS observations, designed to 
satisfy the needs of a variety of users (http://modis-snow-ice.gsfc.nasa.gov/). These 
include a Level-2 swath product; daily and 8-day composite Level-3 tile products which 
are mapped onto a sinusoidal projection and available in 10 degree lat/lon tiles; as well as 
daily, 8-day composite, and monthly Level-3 products available in the Climate-Modeling 
Grid (a latitude-longitude grid) at 0.05 degree resolution (Hall et al. 2005, Hall et al. 
2002, Riggs et al. 2005, Riggs et al. 2006). An 8-day composite is considered useful 
because in many regions, particularly at high latitudes, persistent cloudiness limits the 
number of days available for surface observations (see results section). The Climate- 
Modeling Grid was developed to be useful for the evaluation of climate models and for 
studies at large spatial scales. Other features of the MODIS snow product suite include 
daily snow albedo (Klein and Stroeve 2002) and fractional snow cover (Salomonson and 
Appel 2004). In addition, a new cloud-gap filled product provides information on cloud 
persistence, and uses observations from prior days to map snow (Hall et al. 2010). Figure 
2 shows examples of MODIS snow cover maps in swath format (MOD10_L2) following 
a major snowstorm in the northeastern United States in December 2010. The analysis 
presented in sections 4 and 5 of this article uses the MOD10C1 daily Level 3 global .05 
degree daily Climate-Modeling Grid product with a spatial resolution of ~5km. 
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1 Validation of the suite of MODIS snow cover products has been undertaken by many 

2 authors as described in Hall and Riggs (2007). These products have also been used 

3 extensively in modeling efforts, both at the regional and hemispheric scales. A 

4 bibliography of papers utilizing the MODIS snow cover products for both validation and 

5 modeling may be found at: http://modis-snow-ice.gsfc. nasa.gov/?c=publications. 

6 

7 3.2 passive microwave based products 

8 

9 Historical passive microwave measurements are available from the Scanning 

10 Multichannel Microwave Radiometer (SMMR) instrument (1978 through 1987), and the 

11 Special Sensor Microwave / Imager (SSM/I) instrument (1987 through present) although 

12 some compatibility issues between the two products exist, due to slight differences in the 

13 frequency bands measured, overpass time, swath width, native footprint resolution, and 

14 coverage issues related to SMMR being powered down every other day. (Armstrong and 

15 Brodzik 2001, Brodzik et al. 2007, Derksen and Walker 2003). AMSR-E, available from 

16 2002 to October 2011, provides a suite of measurements to make it spectrally compatible 

17 with both SMMR and SSM/I at higher spatial resolution (Derksen et al. 2005b, Kelly et 

18 al. 2004). Due to the inherent difficulties and regional variations in the interpretation of 

19 passive microwave signals, the production of a data set that is consistently accurate 

20 across all Northern Hemisphere regions requires either (1) a physical approach, which 

21 includes robust representations of snow pack processes and their parameterization in 

22 retrieval schemes (Pulliainen and Hallikainen 2001), or (2) a regional approach, which 

23 includes regionally-tuned algorithms (Foster et al. 1997) that statistically represent 
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1 regional snow pack processes but which are not applicable in different snow 

2 accumulation regimes. The physical approach is very challenging yet has the potential of 

3 being widely applicable as our knowledge of, and ability to model, regional snow pack 

4 processes improves, while the statistical approach is applicable only in the few regions 

5 for which retrieval schemes have been calibrated. 

6 

7 The global AMSR-E SWE product suite (Tedesco et al. 2011 updated) consists of daily, 

8 pentad (five-day) maximum and monthly average SWE estimates that together comprise 

9 the only NASA satellite-based SWE product available to the scientific community. As an 

10 example, Figure 3 shows the snow depth obtained from the AMSR-E product for January 

11 30 th , 2005. 

12 

13 The AMSR-E SWE operational algorithm takes advantage of several AMSR-E channels that 

14 are unavailable from SSM/I and SMMR. Snow depth is derived as a combination of 

15 microwave brightness temperature differences at different frequencies, weighted by 

16 coefficients derived from the difference between vertical and horizontal polarizations. These 

17 coefficients replaced a previously used static coefficient to attempt to capture the spatio- 

18 temporal variability of parameters such as grain size (Kelly 2009, Kelly et al. 2003, Tedesco 

19 and Narvekar 2010). The algorithm uses a spatially varying but temporally static map of 

20 snow density. 

21 

22 Environment Canada (http://ccin.ca/cms/en/socc/snow/swe/currentSnow.aspx) also 

23 produces a regional passive microwave SWE product for central Canada, including the 

24 Prairies and part of the boreal forest back to 1978. Until December 1999, this product 
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1 relied on a single algorithm that was calibrated for the prairies region (Goodison and 

2 Walker 1995). Since that time algorithms that correct for the effects of different forest 

3 types (Derksen 2008, Goita et al. 2003) and the sub-Arctic tundra (Derksen et al. 2010) 

4 have been developed . 

5 

6 3.3 Combined Products 

7 One promising avenue, and an area where great efforts are currently being made, is to 

8 refine our abilities to combine ground observations and models with space borne 

9 remotely sensed data. Tait et al. (2000) provide a helpful review, and describe a prototype 

10 of a fully automated product that includes station observations as well as both visible and 

1 1 microwave retrievals. Here we briefly review some products that include combinations of 

12 satellite, station observations, and models. While not exhaustive, it provides examples of 

13 the variety of integrative products and methods that have been produced. 

14 

15 3.3.1 The Canadian Meteorological Center snow product (station observations and 

16 modeling ) 

17 Snow depth from the Canadian Meteorological Center Daily Snow Depth Analysis 

18 includes a hybrid modeling / observational approach based on optimal-interpolation of 

19 daily snow depth observations from hundreds of stations globally, with snow density 

20 estimated from a simple snow pack model (Brasnett 1999). This model output is 

21 considered most dependable over regions with significant station coverage, which is 

22 generally south of 55° North, where model results are well constrained by observations. 

23 Over most of the Arctic, in contrast, where there are few observations, the analysis is 
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based mostly on model results, and is skewed towards snow depth observations at coastal 
locations with observing sites at open areas near airports. Snow at these sites tends to be 
shallower and to melt out earlier than snow in surrounding terrain. Nevertheless, this 
analysis is considered to be a reasonable estimate of snow depth over data-poor Arctic 
regions, and has been used in a number of studies (Brown and Mote 2009). Here we use 
CMC modeled snow depths for comparison with AMSR-E snow depths. 

3.3.2 GlobSnow (satellite, station, and model) 

In 2008 the European Space Agency embarked on an effort to develop a long term snow 
cover data set called GlobSnow with sufficient homogeneity to be acceptable for climate 
change analysis. The GlobSnow product currently includes global gridded information on 
snow extent and SWE across the Northern Hemisphere (excluding mountainous regions) 
(Pulliainen 2010). The SWE product is based on the method of Pulliainen (2006). By 
incorporating station observations and snow pack modeling into passive microwave 
retrieval algorithms, the goal is to provide an accurate product useful for analyses at 
many different spatial scales, and for near-real time as well as climatological studies. The 
snow extent product is created using European Space Agency satellite visible and 
infrared observations (ERS-2 ATSR-2 and Envisat AATSR) based on the method of 
Metsamaki (2005), and will likely be available at a variety of spatial resolutions. 
GlobSnow is currently available ( http://www.globsnow.info/) but is new, so there is little 
peer-reviewed literature on it at the time of this writing (Takala et al. 2011). 
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1 3.3.3 Other combined products 

2 A variety of combined products have been produced globally, regionally, or for specific 

3 purposes. One widely used combined product is NOAA’s National Operational 

4 Hydrologic Remote Sensing Center Snow Data Assimilation System, which operationally 

5 incorporates input from snow models, station reports, and airborne sensors to estimate 

6 daily SWE at 1 km resolution across the continental US (Carroll et al. 2001). The product 

7 by Brown et al. (2003), which employs the operational snow depth routine of the 

8 Canadian Meterological Center model (Brasnett 1999), has been used for evaluation of 

9 climate models (Frei et al. 2005). Foster et al. (2011) recently produced a global product 

10 blending visible and infrared, passive microwave, and active microwave scatterometer 

11 data, with the intention of incorporating the most reliable aspects of each product. 

12 Derksen et al. (2004) produced a product going back to the early 20 th century for the 

13 North American Prairies and Great Plains based on passive microwave and station 

14 observations. Biancamaria et al. (2011) estimated Northern Hemisphere fields of SWE 

15 based on passive microwave combined with a dynamic snow grain model. Grundstein et 

16 al. (2002) developed a research-oriented SWE climatology for the Great Plains of the 

17 United States by combining station observations with the 1 -dimensional snow pack 

18 model SNTHERM (Jordan 1991). A research-oriented product based on spatial 

19 interpolation of in situ depth measurements over North America (Dyer and Mote 2006) 

20 has been used for process studies (Ge and Gong 2008). The QuickSCAT active 

21 microwave scatterometer has been used to estimate the timing of snow melt across 

22 Greenland (Nghiem et al. 2001) and across Arctic lands (Wang et al. 2008). 

23 
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1 4. Methodology to compare and contrast products 

2 In this section we describe the methodology that we use to demonstrate the regions over 

3 which the products typically differ. This analysis is not meant to provide insight into new 

4 remote sensing techniques, but rather to demonstrate the spatial extents and magnitudes 

5 of the differences between products during different seasons. The methodology employed 

6 here is designed to achieve two goals: (1) to identify regions across the Northern 

7 Hemisphere where there is agreement/disagreement between the three main products 

8 discussed here during clear days; and (2) to provide an indication of the spatial 

9 distribution of uncertainty in the AMSR-E snow depth estimates, as determined by 

10 comparison to the CMC product. In this report we show results for three months: October 

11 (a month of rapid average increase of snow area), January (the month of largest average 

12 snow area), and April (a month of rapid average decrease of snow area). For our analysis, 

13 AMSR-E SWE values are converted to depth. This is done using a fixed density mask, 

14 which is also used as part of the standard product algorithm to estimate SWE values. We 

15 reverse the process in order to convert SWE values to depth. The reprojection methods, 

16 and the methods for each goal, are discussed in more detail in sections 4.1 through 4.3. 

17 

18 Many of our methods for goal 1 closely follow Frei and Lee (2010), and the reader is 

19 referred to that article for more details and justification of the methods. Note that, without 

20 independent verification, agreement between products does not guarantee that they are 

21 correct; and, that if two of the three products agree, it does not guarantee that the third 

22 product is incorrect. 

23 
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1 4.1. Reprojection procedure 

2 IMS and MODIOCI data sets were reprojected to the EASE-Grid 25 km projection 

3 (Brodzik and Knowles 2002) (AMSR-E is already in this projection). Each EASE-Grid 

4 cell value was calculated as a binary (i.e. snow or no-snow) value. Because reprojection 

5 can introduce spurious errors at the grid-cell scale, and these errors are likely to be 

6 exacerbated in areas of variable terrain, we show no results for EASE-Grid cells within 

7 which the GTOPO 30 DEM elevation field has a standard deviation >100 meters. We 

8 also avoid drawing conclusions from individual grid points, but rather focus on results 

9 over large regions with relatively little topographic variation. The reprojection, binary 

10 snow value calculation, and terrain masking were performed according to the method of 

1 1 Frei and Lee (2010). 

12 

13 4.2. Agreement/disagreement between IMS, MODIS, and AMSR-E snow extent 

14 Since both IMS and MODIOCI provide binary values indicating either the presence or 

15 absence of snow (the standard MODIS products also provide fractional snow cover) but 

16 not snow depth, AMSR-E snow depths were converted to a binary value to facilitate this 

17 comparison. All AMSR-E depth estimates below 5 cm are considered snow free as that is 

18 the depth value assigned to shallow snow. 

19 

20 All snow extent analyses include, at each grid cell, only days with “clear” skies, and only 

21 days for which all three products have non-missing data. We use the MODIOCI cloud 

22 mask to identify EASE-Grid cells that are mostly clear. Because MODIOCI. 05 degree 

23 cells are higher resolution than the EASE-Grid and includes fractional cloud cover, they 
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1 can be used to estimate fractional cloud cover within each grid cell of our analysis. And, 

2 because the MOD IOC 1 cloud mask is considered conservative (in the sense that cloud- 

3 covered scenes are unlikely to be designated as “clear”) (Riggs and Hall 2002), we feel 

4 confident that information from cloudy days is not being retained for analysis. This is 

5 achieved by retaining for analysis, for each EASE-Grid cell, only days with >80% of 

6 MOD 10C 1 cells that are <20% cloud covered, for which no product is missing data. Frei 

7 and Lee (2010) present the rationale for this method and explain how the results are 

8 insensitive to reasonable values of these parameters. For each grid point on each day, 

9 either all three products agree (i.e. either snow or no-snow), or one product differs from 

10 the other two. The figures summarizing our results show, for each month, where each 

1 1 product disagrees with the other two products. 

12 

13 4.3. Comparison ofAMSR-E and Canadian Meteorological Centre (CMC) snow products 

14 For passive microwave data and the CMC model, no cloud mask is invoked, so we retain 

15 for analysis all available dates. While passive microwave data are not limited by most 

16 clouds, clouds with high liquid water content can affect the comparison between 

17 spaceborne- and ground-based SWE estimates (Wang and Tedesco 2007); this issue is 

18 ignored here in order to increase the sample size. 

19 

20 The comparison of AMSR-E to the CMC snow product is done by comparing 

21 climatological maps (2003 - 2010). For each month, three panels are shown containing 

22 maps of AMSR-E snow depth, CMC snow depth, and the difference between the two 
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1 products (we calculate the difference as CMC minus AMSR-E, so that a negative 

2 difference indicates that AMSR-E overestimates snow depth with respect to CMC). 

3 

4 5. Results 

5 In this section we show the results of our analysis, the purpose of which is to demonstrate 

6 the spatial patterns of disagreement between the data sets. We also discuss possible 

7 reasons for disagreements. In some cases these reasons may be speculative. 

8 

9 5.1 Number of days per month available for analysis 

10 Before discussing disagreements between the products, we first show maps of the number 

11 of days per month available for comparison (figure 4) which demonstrate the problem 

12 presented by clouds. During October and January (figures 4a,b) most Arctic land surfaces 

13 are colored green or dark blue, which indicates that on average less than three (green) or 

14 three to six (dark blue) days per month are available for analysis. (In January one also 

15 sees the "ring" around the Arctic with no data associated with no solar illumination.) 

16 During spring, which tends to be less cloudy over most regions (figure 4c), one can find 

17 large portions of the Arctic with either six to nine or nine to twelve days per month 

1 8 available for comparison. 

19 

20 The vast majority of the unavailable days are caused by clouds, not by data that is 

21 missing for some other reason. Any satellite product based on visible and infrared band 

22 radiances will lack information from the surface under clouds. While passive microwave 

23 based products can provide information under most types of clouds, they are currently 
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1 unreliable under a number of circumstances (discussed in the next section). Considering 

2 the importance of having daily real-time information about the surface to specify 

3 boundary conditions in weather prediction models, as well to track climatological 

4 changes in snow extent, IMS, or an equivalent product that provides information for all 

5 days regardless of cloud conditions, is a necessity. 

6 

7 5.2 Disagreement between AMSR-E and the other two products 

8 Figure 5 includes, for each month, a map showing where AMSR-E identifies snow to the 

9 exclusion of the both MOD10C1 and IMS (figures 5a,c,e) and a map showing where 

10 AMSR-E finds no snow when the other two products identify snow (figures 5b,e,f). The 

11 most prominent feature is the red colored plateau region of central Asia seen in all maps 

12 down the left hand column (figures 5a,c,e). This indicates that during all months over this 

13 region AMSR-E identifies snow more often than the other two products. While we do 

14 not, in general, assume that a product is wrong because it disagrees with the other two 

15 products, in this region we know from other studies that AMSR-E observations are biased 

16 due to problems in passive microwave snow detection at higher elevations associated 

17 with atmospheric influences on passive microwave retrievals (Wang and Tedesco 2007). 

18 Since the atmosphere over the high elevation plateaus is much thinner, the algorithms 

19 calibrated globally at lower elevations require correction (Savoie et al. 2009). 

20 

21 Panels on the right side of figure 5 show that during each month there are regions where 

22 AMSR-E identifies snow less frequently than MOD10C1 or IMS (figure 5b,d,f). The 

23 regions shown on these panels coincide with boundary of the continental snow cover 
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during each month (see the Rutgers University Global Snow Lab web site for 
climatological maps of monthly snow cover based on IMS: 
http://climate.rutgers.edu/snowcover/index.php) . Regions near the boundary tend to have 
patchy, shallow snow packs. During spring (figure 5f) the disagreement across well 
defined ablation bands at the southern boundary of the continental snow pack is also 
likely due to significant areas of melting snow with liquid water in the snowpack . 

5.3 Disagreement between IMS and the other two products 

Figure 6 demonstrates that the most prominent circumstance under which IMS disagrees 
with the other two products is during the spring ablation period near the boundary of the 
continental snow pack (figure 6e). This result is in agreement with recent studies (Brown 
et al. 2007, Brown et al. 2010, Frei and Lee 2010) which find that over the last decade or 
so the timing of spring ablation over North America is later, by up to several weeks in the 
central Canadian Arctic, according to IMS in comparison to other observations. The 
reasons for these discrepancies, which are found during the entire spring ablation season 
(April, May, and June; May and June not shown here) over the boreal forest as well as 
the tundra, are not understood, but may be related to geographic factors such as the forest 
type and/or the presence of numerous lakes in this (Derksen et al. 2005a, Rees et al. 
2006), Investigations into the cause of this problem continue. 

5.4 Disagreement between MOD10C1 and the other two products 

The most interesting example of MOD10C1 disagreeing with IMS and AMSR-E is found 
during autumn over Eurasia. During October over a broad, seemingly incoherent region 
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1 of Eurasia, predominantly over Scandinavia and northern Europe, MOD10C1 often 

2 identifies snow when the other two products do not (figure 7a). However, this region is 

3 not as incoherent as it may seem, as it corresponds closely to the boreal evergreen 

4 needle leaf forest as defined by analysis of MODIS reflectances (Friedl et al. 2002). 

5 During November (not shown) we find a similar pattern, except the differences are more 

6 extreme and concentrated more over Scandinavia. The eastern Eurasian region over 

7 which MOD10C1 often fails to identify snow when IMS and AMSR-E see snow (figure 

8 7b) corresponds closely to the region of deciduous needleleaf forest. It seems that over 

9 one type of forest MODIS sees snow more often, while over a different type of forest 

10 MODIS sees snow less often. While the difficulties of remotely sensing snow under 

11 forest canopies have been widely reported, the patterns reported here have not been 

12 examined in the literature. 

13 

14 5.5 Comparison of AMSR-E to the CMC snow product 

15 Maximum October snow depth values over the Northern Hemisphere are - 20 - 30 cm 

16 (figure 8). The AMSR-E product suggests more snow in Siberia than the CMC product. 

17 AMSR-E overestimation with respect to CMC over Siberia increases as the snow season 

18 progresses. In January, snow depth differences between the two products increase to - 

19 30-40 cm (figure 9). In April, the area over which AMSR-E overestimates snow depth 

20 increases even further with respect to January. In contrast, over other regions AMSR-E 

21 tends to underestimate snow depth with respect to the CMC product, but these areas do 

22 not appear to expand as the snow season progresses. These include the Tibetan plateau 

23 and along the north-east coast of North America (figure 10). 
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2 Histograms of the snow depth differences between the two products are shown in figure 

3 11. Overall, AMSR-E tends to overestimate the values provided by CMC. While the 

4 variance of the errors can be seen in the histogram plots, perhaps a more informative 

5 number would be the coefficient of variation (Cv). Cv, defined as the absolute value of 

6 the standard deviation of the differences divided by the mean CMC snow depth, provides 

7 an indication of how large the differences are in comparison to the snow depth. For 

8 example, a value of Cv=l means that the errors are of the same magnitude as the mean 

9 depth; Cv=0.1 means that the errors or an order of magnitude less than the mean snow 

10 depth. Cv values were calculated for each month (table 1). Cv values are highest in 

11 October, when depths are small; lowest in January; and increase again in April. As a 

12 snow pack ages, even under cold conditions without additional precipitation, 

13 metamorphic processes lead to grain size variations (such as depth hoar formation) that 

14 tend to introduce errors in the passive microwave product. Furthermore, as temperatures 

15 fluctuate and additional precipitation events add fresh snow, snow packs can develop a 

16 series of well defined layers of different grain sizes that confound passive microwave 

17 based estimates of depth and SWE. Ice layers, which can develop as a result of melt- 

18 freeze and rain-freeze events, introduce additional scattering and therefore additional 

19 uncertainty. Such complications, combined with the impact of vegetation, especially 

20 vegetation that can change seasonally, can introduce a growing error in passive 

21 microwave retrievals as the snow season progresses. 

22 
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1 Improved confidence in our abilities to estimate snow mass from satellites would support 

2 efforts to monitor the fresh water flux into the Arctic Ocean. An order of magnitude 

3 estimate suggests that the volume of water in the snow pack can play a significant role in 

4 the total annual river runoff into the Arctic Ocean of 4300-4800 km 3 yr’ 1 (Arnell 2005, 

5 Miller and Russell 2000). Our AMSR-E based (highly uncertain) estimate of the mean 

6 snow mass over land surfaces during March (the month of maximum snow mass) north of 

7 60 N is -1600 km 3 . Frei et al. (2005) based on the analysis of Brown et al. (2003) 

8 estimated the observed mean snow volume over North America during March to be 

9 -1500 km which was equal to the median value estimated by a group of 18 climate 

10 models. This compares to a recent passive microwave-based estimate of -1400 km and 

11 -2300 km for mean North American and Eurasian snow volumes, respectively 

12 (Biancamaria et al. 2011). The errors associated with most of these estimates are 

13 currently unknown, but they indicate that the snowpack provides a significant fraction of 

14 the total river runoff to the Arctic. 

15 

16 

17 6. Discussion and Conclusions 

18 

19 For most of the snow season and most regions there is large-scale agreement amongst the 

20 products in identifying the location of snow covered surfaces (i.e. snow extent, regardless 

21 of snow depth) during clear sky conditions. One exception to this is over central Asia. It 

22 is known that passive microwave products identify snow on the Tibetan Plateau and 

23 surrounding mountains when visible and infrared products do not (Armstrong and 

24 Brodzik 2001, Basist et al. 1996). Because passive microwave retrieval algorithms are 
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1 calibrated at lower elevations, at these high elevations the reduced atmosphere between 

2 the surface and radiometer can result in retrieval errors (Savoie et al. 2009, Wang and 

3 Tedesco 2007). The second exception occurs where snow is ephemeral, patchy, or wet. In 

4 such regions the attenuation of the passive microwave signal, upon which snow detection 

5 is based, is compensated for by emission from the surface or from liquid water in the 

6 snow pack (Matzler 1994, Ulaby and Stiles 1980). Despite these difficulties, all estimates 

7 (discussed in the preceding section) indicate that the snowpack is the source of a 

8 significant portion of runoff into the Arctic basin. 

9 

10 The disagreements in snow extent during April are greater than during October or 

11 January in terms of the percentage of available days during which one product differs 

12 from the other two. This is in agreement with Brown et al. (2010), Frei and Lee (2010), 

13 and Brown et al. (2007), who find differences between sensors during spring over North 

14 American regions experiencing ablation, and indicates that the wet snow during ablation 

15 is perhaps more of a hindrance to the identification of snow from satellites than some of 

16 the other confounding factors. However, during fall and winter the evaluation is 

17 hampered by data availability problems associated with cloudiness and solar illumination 

18 issues. 

19 

20 Our analysis also demonstrates that snow depths estimated by the Canadian 

21 Meteorological Centre product and by the AMSR-E algorithm can differ substantially. 

22 Although there are no absolute surface reference observations in most regions to 

23 determine which (if either) product is correct, we know from experience as well as theory 
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1 that the passive microwave depth and SWE algorithms are inaccurate under certain 

2 conditions (Tedesco and Narvekar 2010). Sources of error include: surface heterogeneity 

3 within a passive microwave footprint; temporal and spatial variability in grain size and 

4 snow density; obscuration of snow by forests; masking of the passive microwave signal 

5 by liquid water in the snow pack; and effects of atmospheric attenuation. The persistent 

6 underestimation by AMSR-E with respect to CMC over some regions can be partially 

7 explained by considering that snow depth over many of those areas is above the 

8 ‘saturation’ depth to which the passive microwave algorithm is sensitive (Derksen 2008, 

9 Markus et al. 2006, Matzler 1994, Schanda et al. 1983); the presence of a high fraction of 

10 lakes over the north east of North America is also believed to be a source of error 

11 (Derksen et al. 2005a, Rees et al. 2006). 

12 

13 Another example is the overestimation of snow depth by AMSR-E over northern Siberia, 

14 which can be attributed to the limitation of the current AMSR-E algorithm to account for 

15 the large grains that typically develop in snow packs in this (and some other) regions 

16 (Clifford 2010). Over regions that develop and maintain a snow pack early in the season, 

17 the snow tends to insulate the ground keeping it warm even as air temperatures fall, 

18 resulting in a strong vertical temperature gradient in the snow pack. This temperature 

19 gradient causes vertical energy and vapor fluxes within the snow pack, the net effect of 

20 which is a layer of depth hoar at the bottom of the snow pack (Jordan et al. 2008). The 

21 large crystal sizes of depth hoar (~5mm) cause increased scattering of microwave 

22 radiation resulting in an overestimation of the snow pack by the passive microwave 

23 algorithms. 
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1 


2 Opportunities remain for the development of improved snow products. For example, 

3 improvements can be made with regard to the retrieval of snow amount from passive 

4 microwave sensors (Tedesco et al. 2004) under forested terrain (Derksen 2008), the 

5 refinement of snow extent estimates from visible and infrared sensors (Parajka and 

6 Bloschl 2008), and the estimation of sub-grid scale information. Tedesco and Miller 

7 (2007b) explore the relative merits of combining active and passive microwave retrievals, 

8 using a MODIS snow product as their reference “truth.” A number of researchers are 

9 investigating the potential for finer scale information on snow extent, amount, fractional 

10 snow cover (Derksen et al. 2005b, Salomonson and Appel 2004, Salomonson and Appel 

11 2006), snow melt (Wang et al. 2008), as well as on snow pack properties (Kinar and 

12 Pomeroy 2007, Nolin and Dozier 2000, Painter and Dozier 2004, Painter et al. 2003, 

13 Rango et al. 2000, Schmugge et al. 2002). Improvements in remotely sensed products 

14 that do not rely on the assimilation of data or model results will come as a consequence 

15 of improved understanding of the interaction between electromagnetic and geophysical 

16 parameters at large spatial scales. In this context, a new operational algorithm based on 

17 the inversion of an electromagnetic model, artificial neural networks and snow 

18 climatology currently under evaluation may be capable of accounting for some of these 

19 limitations. 

20 

21 One currently active area of research is the development of combined products, which 

22 include in situ observations and/or modeling results as well as remotely sensed 

23 information. One can identify advantages and disadvantages to both combined and Stand- 
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1 alone remotely sensed products. While stand-alone remotely sensed products contain 

2 inherent drawbacks as discussed here, at any time, either in situ or remotely sensed data 

3 streams can fail, rendering combined as well as stand-alone products vulnerable to 

4 missing information. This is most critical for real- or near-real time operational products, 

5 on which weather forecast models or time- sensitive decisions rely. 

6 

7 Remote sensing of snow continues to contribute to our understanding of Earth system 

8 processes. MODIS snow products are valuable because they can provide high resolution 

9 snow estimates under cloud-free conditions using a quantifiable algorithm. However, for 

10 climatological as well as operational purposes, humans can integrate and filter data from 

11 multiple sources and satellite images in ways that fully automated methods are (at least 

12 currently) unable to, and provide information for the entire land surface of the globe, 

13 regardless of the presence of clouds. Thus, continuation of IMS, with its long record of 

14 snow extent, is a priority. Considering the difficulties in determining SWE on a global 

15 scale from stand-alone remote sensing products, it seems likely that combining multiple 

16 sensors with station observations and/or models, such as in the GlobSnow product, will 

17 provide the best estimates of SWE. 

18 

19 
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1 Table 1. Mean snow depth from CMC product; standard deviation of the differences 


2 between the CMC and AMSR-E snow depths; and the coefficient of variation. All values 

3 are averages of grid points across all Northern Hemisphere land areas north of 30 N 

4 excluding the Greenland ice sheet. 

5 



p (mean CMC 
snow depth) 

(cm) 

a (standard 

deviation of 

difference) 

(cm) 

Cv abs(a/p) 

(coefficient of 

variation) 

(unitless) 

October 

5.0 

4.15 

0.83 

January 

22.5 

7.18 

0.32 

April 

22.9 

12.33 

0.54 
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